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We have performed surface tension measurements on carbopol gels of different concentrations and 
yield stresses. Our setup, based on the force exerted by a capillary bridge on two parallel plates, 
allows to measure an effective surface tension of the complex fluid and to investigate the influence 
of flow history. More precisely the effective surface tension measured after stretching the bridge 
is always higher than after compressing it. The difference between the two values is due to the 
existence of a yield stress in the fluid. The experimental observations are successfully reproduced 
with a simple elasto-plastic model. The shape of successive stretching-compression cycles can be 
described by taking into account the yield stress and the elasticity of the gel. We show that the 
surface tension ynv of yield stress fluids is the mean of the effective surface tension values only if the 
elastic modulus is high compared to the yield stress. This work highlights that thermodynamical 
quantities measurements are challenged by the fluid out-of-equilibrium state implied by jamming, 
even at small scales where the shape of the bridge is driven by surface energy. Therefore setups 
allowing deformation in opposite directions are relevant for measurements on yield stress fluids. 


I. INTRODUCTION 

Yield stress fluids are widespread materials in every¬ 
day life, food industry, cosmetics, building industry, oil 
industry and many other fields. They are of a great in¬ 
terest because they have the property to flow only when 
the applied stress is greater than a critical stress called 
yield stress [T]. They include emulsions, suspensions, 
gels, granular pastes and foams. 

Their bulk properties have been studied extensively 
since the work of Herschel and Bulkley in 1926 [ 2 ] and 
are now well characterized [3]. Besides, capillarity and 
wetting are well known for simple fluids. Recently a lot of 
work has also been done on surface tension and wetting of 
soft solids mg and on the competition between capillary 
forces and elasticity of the substrate BIS]. But until 
now few studies have focused specifically on the surface 
tension of yield stress fluids [S]. 

However surface tension and wetting properties of yield 
stress fluids are of a great importance for capillary imbi¬ 
bition, coating, surface instabilities and adhesion, among 
other applications. 

Here we explore the competition between surface ten¬ 
sion, which is an equilibrium property to be measured, 
and yield stress effects that often keep the system out of 
thermodynamical equilibrium due to a dynamical arrest 
of flow. This situation can be compared to contact angle 
hysteresis: the contact angle is always smaller or greater 
than the equilibrium (Young) value, depending on the 
history of the contact line m- 

Geraud et al. studied this competition in capillary rise 
experiments m-, this method allowed to measure the 
surface tension and the yield stress of the fluid at the 
same time. Yet, with their setup, a large amount of liq¬ 
uid is needed, the contact angle must be measured in 
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another experiment, the results are extremely sentitive 
to the least defect on the inner surface of the capillar¬ 
ies and only fluids of low yield stress (< 20 Pa) can be 
characterized this way. 

The method presented here allows to get rid of these 
difficulties. Moreover both extension and compression of 
the system can be imposed, which highlights the effect of 
the flow history on the effective surface tension measured. 

The article is built as follows. In the first part, we 
present the fluids on which we performed measurements 
and the experimental setup. Then we describe the exper¬ 
imental results and in the next part we compare them to 
an elastoplastic model using only few ingredients. Finally 
we discuss the agreement between the experiments, the 
model and results from other works. 


II. MATERIALS AND METHODS 
A. Fluids 

The simple fluids used here are deionized water (18 
Mfl) and silicon oil (47V100 from Roth). The complex 
fluids are carbopol gels of different concentrations. The 
raw polymer (powder) is ETD 2050 from Lubrizol. The 
gel is prepared as follows: a small amount of polymer is 
weighted and slowly dissolved in deionized water heated 
at 50°C and stirred. The hot solution is stirred for 30 
minutes, then it is let to cool down to room tempera¬ 
ture. Evaporation is hindered by covering the container 
with Parafllm. Sodium hydroxy de (lOM) is added to the 
solution until its pH is raised to 7 ± 0.5, which causes 
the polymer chains to charge negatively. The charged 
chains thus repel each other, the polymer blobs swell and 
jam, and the solution becomes a gel. Finally the gel is 
either stirred gently by hand or stirred for 24 hours at 
2100 rpm with a mechanic stirrer. It was indeed shown 
in other works m that stirring changes the rheology of 
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carbopol. Our carbopol concentrations range from 0.25% 
(in weight) to 2%. Hand stirred (respectively machine 
stirred) carbopol is denoted HS (respectively MS) in the 
following. 


B. Rheology 

Carbopol gels are generally considered as model, non 
thixotropic, yield stress fluids. As long as slip m , tran¬ 
sient shear banding [14] and confinement [15] are avoided, 
their flow curve, relating the shear stress a to the shear 
rate 7 , is well fitted with a Herschel-Bulkley (HB) law: 

(7 = 0 if tr < (Ty 

1 CT = cry -|- K'^ " if tr > (Ty 

with CTy the yield stress, K the consistency and n the HB 
exponent. 

Our rheometer is an Anton Paar MCR 301 equipped 
with a rough cone and plate geometry of angle 4°. The 
flow curve is obtained with decreasing steps of constant 
shear rate, ranging from 100 s“^ to 0.01 s“^ (10 points 
per decade). The duration of each step is set automati¬ 
cally by the rheometer (between 15 s and 30 s per step) 
and the measurement is made when the steady state is 
reached. The elastic modulus G' is measured by oscilla¬ 
tory shear deformation of 1 % with increasing, then de¬ 
creasing frequencies, ranging from 0.1 Hz to 50 Hz. 

The yield stress of our carbopol samples ranges from 
CTy = 0.3 Pa to CTy = 38 Pa and their elastic modulus at 
0.1 Hz from G' = 1.5 Pa to G' = 155 Pa. Typically n is 
always between 0.5 and 0.6 and K ranges from 0.75 Pa.s" 
to 13 Pa.s". 


C. Bridge tensiometer 

1. Setup. 

The home made bridge tensiometer [16] (figure al¬ 
lows to measure the surface tension of fluids. It consists 
in two horizontal glass plates, between which a droplet of 
the liquid of interest (a few microliters) is deposited. The 
liquid forms a capillary bridge between the two plates. 
The bottom plate is attached to a micromanipulator, so 
that its position can be adjusted by the operator. The 
force applied by the bridge on the top plate is recorded 
through a flexible cantilever equipped with an electro¬ 
magnetic deflection sensor. A high resolution camera 
(Pixelink PL-A 686 M) coupled to an optical magnifier al¬ 
lows to take pictures of the bridge. An example of picture 
is shown in figure 

Both top and bottom plates must be perfectly cleaned 
to avoid line pinning which could deform the axisymmet- 
rical bridge, and to avoid polluting the fluid with dust 
or surfactants. Before each series of measurements, the 
bottom plate is always thoroughly cleaned in a plasma 


deflection sensor 



Figure 1. Drawing of the bridge tensiometer setup. Inset: 
Example of picture of a carbopol bridge. The white stain in 
the middle is a deformed image of the fiat LED light situated 
in the back of the setup. 


cleaner. The small top plate is dipped in piranha solution 
(1 part of hydrogen peroxyde in 2 parts of concentrated 
sulfuric acid) and rinsed with deionized water. 



Figure 2. Definition of the main geometrical parameters. 


2. Force balance. 

At equilibrium the force on the cantilever (A) and the 
geometry of the bridge are directly linked via the surface 
tension of the fluid {'yLv)- More precisely the force mea¬ 
sured by the cantilever is the sum of the pressure force 
at the liquid-plate interface and of the capillary force at 
the perimeter of this interface HZ]: 

F = —ttRq X Ap + 2 ttRq sin 00 x iLV 

where Rq and 9q are the radius and the contact angle 
defined on figure assuming cylindrical symmetry, and 
Ap is the pressure difference between the fluid and the 
atmosphere. 

It must be noticed that the same force balance can be 
done at each height z of the bridge, and especially at the 
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neck (zat) where sui0{zn) = 1- However to account for 
gravity, it is necessary to add the weight of fluid above 
Zn (denoted W) to the force balance: 

F = W — ttR]^ X Ap + 27ri?Ar x 

Rn being the radius of the bridge at the neck. 

Finally Laplace’s law allows to replace the pressure dif¬ 
ference Ap with X Cat, Cat being the total curvature 
of the surface at the neck: 

F-W = -fLvx - t^RnCn) = ILV X L (1) 

3. Measurement protocol. 

To form the bridge, a droplet is deposited on the bot¬ 
tom plate which is then moved upwards until contact of 
the liquid with the top plate. Generally the liquid spreads 
on the whole upper plate and the two plates are stuck to¬ 
gether, so the bridge must always be stretched before the 
beginning of the measurement. During the experiment, 
the bridge is stretched or compressed by changing the 
position of the bottom plate and then let to equilibrate. 
Because of evaporation the force is never completely con¬ 
stant, but the force value and the picture are saved when 
the force evolution is sufficiently slow (about 1 /iN per 
second) compared to the total force step (of the order of 
100 fiN in a few seconds). A typical example of force step 
is shown in figure [T6| of Appendix A. 

For each aspect ratio of the bridge, the geometric pa¬ 
rameter L (defined in equation is computed from the 
picture. The outline of the bridge profile is extracted 
from the image and fitted with a high-order ( 11 ) poly- 
nom. It is necessary to get a rather smooth profile be¬ 
cause it is derivated twice to compute the curvature, but 
the polynom must nevertheless follow the real profile as 
faithfully as possible. The curvature of the surface is 
computed as: 

^ yRjz) _ R''{z) 

^ ’ (1 + i?'(z) 2 )l /2 (1 + i?'(z)2)3/2 

The force is obtained from the cantilever deflection af¬ 
ter calibration. The weight W is a small correction to 
the force and it comes from the calculated volume of fluid 
above zat (obtained by integration of the profile) times 
the fluid density. This way, F — W can be plotted as a 
function of the corresponding parameter L. In the case 
of a simple fluid, a proportional relation is expected (see 
relation[^, the slope being the surface tension value 'yLV- 

III. EXPERIMENTAL RESULTS 
A. Simple fluids 

In order to validate the setup, the experiment was first 
performed with simple fluids. The protocol was always 


the same, testing stretching as well as compression to 
check the influence of the dynamics history on the results. 

With pure water and silicon oil, the force-L plot in¬ 
deed shows a proportional relation (see figure and 
the slopes correspond to respective surface tensions of 
(74 ± 1) mN/m and (21 ± 1) mN/m. The expected sur¬ 
face tensions are 73.0 mN/m and 21.0 mN/m (at 18°C). 
The agreement is very good, with precision comparable 
to usual surface tension measurement methods [ISH2I]- 




Figure 3. Force-L plot for deionized water and silicon oil. 
The force F — W is proportional to L and the slope corre¬ 
sponds to the surface tension of the liquid. Red triangles and 
black squares respectively stand for stretching and compres¬ 
sion steps. 

Great care must be taken to avoid hysteresis of the 
contact line. A treatment was applied on the glass plates 
| 22 | to minimize hysteresis before the measurements with 
oil and water. 


B. Carbopol 

Garbopol is a yield stress fluid, i.e. it cannot flow if 
the applied stress is below the yield stress (cry). In par¬ 
ticular it is necessary to make small droplets such as the 
surface tension induced pressure dominates over the yield 
stress. This can be quantified by a dimensionless num¬ 
ber comparing the yield stress energy cry x (r being 
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the characteristic size of the system) and the surface en¬ 
ergy 7iy X For example, in a droplet of gel with 
cry = 20 Pa and — 60 mN/m, the excess pressure 
in the liquid due to capillary forces must be greater than 
20 Pa for the droplet to flow and adopt a spherical shape. 
Therefore we expect that an isolated droplet of this fluid 
of radius smaller than 6 mm (r < 2'yi^y jay) is spherical. 
More generally if the curvature of the surface is greater 
than (Jy j^LVi the shape of the liquid surface should be 
controlled by surface tension. This is the reason why 
only small capillary bridges (radius of the order of the 
millimeter) were studied in this experiment. 


in figure We observe that the points do not all align 
on a single line. The solid red triangles correspond to the 
first series of stretching, starting at the top-right angle of 
the plot. The red line is the linear fit of these points, and 
its slope is denoted j^^y. The black squares correspond to 
the series of compressions. They align on a second line, 
whose slope jly is always smaller than for the stretched 
points. This behavior is reproducible for every sample of 
carbopol, and the greater the yield stress, the wider the 
difference of slopes between the two sets of points. 

120r 




Figure 4. Force-L plot for two different carbopol samples: a) 
0.25% (HS), yield stress cry = 5 Pa ; b) 1% (MS), yield stress 
ay — 19 Pa. The solid (resp. empty) red triangles stand for 
the first (resp. second) series of stretchings, the black squares 
for the series of compressions. The linear fit values are written 
in the figure. 
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Figure 5. Upper (red) and lower (black) slopes of the force-L 
plots plotted as a function of the yield stress. The green line 
is a guide for the eyes, indicating the mean surface tension of 
vanishing yield stress carbopols. The error bars indicate the 
averaged points (see text). The error on all the other points 
is ±5 mN/m. 



(Ty (Pa) 


As for the simple fluids we started with a series of 
stretchings and then a series of compressions. Most of the 
time these were followed by a second series of stretchings. 
It must be clear that compressed bridges correspond to 
large values of L and that a stretching is transposed into 
shift to the left side of the force-L plot. On the contrary 
a shift to the right of the plot is a compression (see figure 

H- 

Typical force-L plots for carbopol 0.25% (HS, cry = 5 
Pa) and carbopol 1% (MS, cry = 19 Pa) are reproduced 


Figure 6 . Difference Aqny = 7 /‘y — 7 /y of the force-L plots 
slopes, as a function of the yield stress of the samples. Each 
point color represents the volume of the droplet. Star-shaped 
points stand for HS carbopol and dots for MS carbopols. The 
line indicates a linear fit whose correlation coefficient (R)^ is 
only 0.74. 

To confirm the influence of the yield stress on the 
apparent surface tension, we performed several experi¬ 
ments, varying cry between 0.3 Pa and 38 Pa. This could 















5 


be achieved by varying either the polymer concentration 
or the stirring. Hand-stirred carbopols have indeed a 
mnch greater yield stress than machine-stirred carbopols 
of same concentration. 

For a few samples the experiment was performed with 
several droplet volumes between 2 /xL and 15 /xL. More¬ 
over for two of them, 10 identical measurements were 
carried out in order to evaluate the dispersion of the ef¬ 
fective surface tension values. The standard deviation of 
the results is of about 5 mN/m for given yield stress and 
volume. 

Figurej^shows the values of the upper and lower slopes 
as a function of the sample yield stress, and each point is 
an average on 1 to 4 droplets of similar volume (within 
1 /xL steps) and yield stress (within 1 Pa steps). It can 
be observed that the upper slope increases with the yield 
stress while the lower slope decreases. For vanishing yield 
stress, they both converge to 63 mN/m. 

Figure is a plot of the slopes difference = 

vs ay, with the same average as before, and 
the droplet volume is represented by the point color. It 
confirms the monotonic dependence of the slopes differ¬ 
ence with the yield stress, and it also shows that greater 
Ajlv often correspond to larger drops, for a given yield 
stress. 

In both figures the star-shaped points stand for HS 
carbopol samples and the other points for MS samples. 
The averaged points are indicated by error bars. 

Finally, as shown on figure the shape of the second 
stretching cycle (empty red triangles) varies from one ex¬ 
periment to another. The second stretching set of points 
joins the first stretching line (red) faster when the elastic 
modulus of the carbopol is higher, for equal yield stresses. 


IV. ELASTOPLASTIC MODEL 

To understand the influence of the different parameters 
in our experiment, we have developed a simple model. 
The goal is to understand the role of the flow history on 
the curves obtained with a yield-stress fluid. 

Because the experiments clearly show an influence of 
both the yield stress and the elasticity of the fluid, we 
consider an elastoplastic fluid: below ay it behaves as an 
elastic solid, and at ay it flows until it reaches a station¬ 
ary state. We neglect the consistency K of the Herschel- 
Bulkley model as the time evolution of the force is not 
investigated here, only the final state. The model mimics 
the experimental protocol and explores the influence of 
the elastic deformation on the stress state of the bridge 
for either stretching or compression and different initial 
conditions. To be able to calculate the stress, we consider 
two limiting simplified geometries, the filament and the 
pancake. 

This model allows to faithfully reproduce the experi¬ 
mental results and thus to explain the observations ex¬ 
posed in part |III| as explained below. 

A drop of viscoplastic liquid with yield stress ay and 




Figure 7. Example of two force-L plots for carbopols of same 
yield stress ay = 7 Pa and different elastic moduli, a) MS 
carbopol, G' = 20 Pa. b) HS carbopol, G' = 45 Pa. 


shear elastic modulus G' is considered. The drop has a 
nearly cylindrical shape with height h and neck radius 
Rn, so that the volume of the drop is H « TTR'fjh. We 
denote 9q = 30° the contact angle, which is roughly the 
contact angle observed in the experiments. The total 
curvature is assumed to be constant along z, and the 
geometric parameter L is approximated by 


L — ‘I'KRjy — TrR'jq ( —- 

Hn 


2 cos 6q 


= ttR 


N ■ 

V 

7rh 


27ri?^ cos 00 
h 

2V cos 00 


For a given volume V, the filament (resp. pancake) 
geometry corresponds to heights h ^ {V/ttYA (resp. 
h ^ (V/tt)^/^). The volume is fixed to V = 10 mm^, as 
often encountered in experiments. This corresponds to 

(V/7r)i/3 ~ 1.5 

mm. As this is the typical experimental 
value of h, the experiments do not correspond to any of 
these limiting geometries (filament or pancake), but to an 
intermediate regime where h ^ R. But as discussed later, 
we show that the results of the model do not qualitatively 
depend on the chosen geometry. This is the reason why 
we choose to explain the model in details in the filament 
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geometry only. The pancake geometry calculations are 
nevertheless presented in Appendix B. 


A. Filament geometry 

In this geometry, usually encountered in capillary thin¬ 
ning or filament-stretching devices |231 12^ . elongational 
deformation and normal stress (and not shear) are as¬ 
sumed to be dominant. The characteristic stress is the 
normal stress a = — arr- We look at small height 

variations Ah. The corresponding step in deformation 
is: 



and the stress before each step is denoted cto- 

In the elastoplastic hypothesis, and taking into account 
the tensorial formulation of the stress tensor |^, the new 
stress after a step is given by the following function: 

( if CTo -I- SG'Ae < —y/Say 

^ \ cTo + SG'Ae if —\/Z(Jy < cto + 3G'Ae < -s/Scry 

\ if ^0 T SG^Ae: > -t-y^fTy 

Finally, the normal elastoplastic force applied on the 
cantilever is evaluated at each step: 

2 ^ 

Fep = (T7ri?jv = 

In the experiments the drop is initially stretched so the 
initial stress is set to -t-'x/Scry. Then successive steps of 
deformation Ah = 0.3 mm are applied to the model drop, 
starting with stretching from h = 1.5 mm to h = 4.5 mm, 
then compressing and finally stretching again. 

For each step, the total traction force, which is the 
sum of the capillary force ^lvL and the elastoplastic 
one Tep, is calculated for = 60 mN/m and cry = 
2 Pa or 5 Pa, using the following approximation of L = 
TTyJVliTTh). Several values of the elastic modulus G' are 
tested: G'jay = 0.5, 2 and 8 . 

On figure plots a) and b) differ only by the yield 
stress value. It is clear that the slopes difference be¬ 
tween the two branches increases with the yield stress 
cry. Here the L-range is small compared to the experi¬ 
ments because elongated shapes correspond to low values 
of L. 

Figure [^represents three typical force-L plots from the 
model, for a given yield stress (try = 5 Pa) and differ¬ 
ent elastic moduli G'. It shows that the elastic modulus 
has a strong influence on the shape of the stretching- 
compression cycle. First the shape of the lower branch 
(black symbols) changes when G'/try increases, causing 
a variation of the y-intercept and of the slope of the black 
linear fit. Then the second stretching branch (empty tri¬ 
angles), as in the experiments, joins the first stretching 
branch (solid triangles) all the faster as G' is large. The 




Figure 8. Results from the model: F as a function of L for 
a filament geometry, with G' joy = 8. a) cry = 2 Pa. b) 
CTy = 5 Pa. 


first stretching branch does not change because the nor¬ 
mal stress is -bv^cry all along and G' plays no role here. 

Note that for values of G' jay of the order of 10 or 
more, the maximal elastoplastic stress (-x/Scry) is reached 
immediately after the direction change. This means that 
for G'/cry > 10 the points fall on two limiting curves 
determined only by the yield stress. These two curves 
are symmetrical with respect to F — W = This 

allows to find the true value of y^y by taking the mean 
of the two limiting slopes. 


B. Pancake geometry 

We also checked the other limit of a flattened drop. 
In this case, the deformations and dissipation are domi¬ 
nantly due to shear along the z direction. Therefore we 
cannot use a homogeneous description but we need to 
describe the stress profile at the wall. 

The details of the calculation in the lubrication ap¬ 
proximation can be found in the appendix B. 

What is observed with this geometry does not differ 
qualitatively from the case of the filament (see figure [l0| : 
stretching (resp. compression) of the capillary bridge is 
associated with an increase (resp. decrease) of the ap- 
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Figure 9. Results from the model: F as & function of L for 
a filament geometry, ay = 5 Pa and G'/ay = 0.5, 2 and 8 
from a) to c). 


parent surface tension. The L-range increases, which is 
consistent with the displacement towards the right of the 
plot with compression. 


V. DISCUSSION 

A. Influence of yield stress and volume 

The model confirms the influence of yield stress on the 
difference of effective surface tensions. For large enough 
elastic moduli G' ay the excess force due to the yield 


Figure 10. Results from the model: F as a function of L for 
a pancake geometry and ay = 10 Pa. a): G'/ay = 0.2. b): 
G'/ay = 2. 


stress can be approximated by x L. In the fil¬ 

ament geometry, assuming that the stress has reached 
its saturation value, this excess force can be estimated 
by 'Jiay X and L « 7ri?jv so the slopes differ¬ 
ence reduces to cx R]\[(Ty. In the pancake geom¬ 

etry, the excess force is about ^ay x R%/h and 
L « 2-kR?^ cos O^/h so the relation Ay^y oc Rn^y still 
holds. 


To refine the interpretation, we rescaled our experi¬ 
mental data with the droplet size. Namely, considering 
the most compressed state (indicated with an asterisk), 
the yield stress was multiplied by the neck radius R/y. 
The effective surface tension difference Ay^y shows to 
be proportional to the resulting quantity. The alignment 


of the data points is better after rescaling (figure 11 
= 0.85) than for the raw data (figure = 0.74) 

and the prefactor is of order 1. 

This evidences that even a static surface tension mea¬ 
surement will depend on the flow history, and this all the 
more as the yield stress is high and the droplet is large. 
The error on the measurement, if it is performed after 
(or during) a flow in always the same direction, will be 
of the order of ay x r with r a dimension of the system. 

The length r must be thoroughly identified. In our 
experiments the bridge radius at the most compressed 
















state R’^ is the characteristic length scale because it cor¬ 
responds to the greatest force difference in a force-L plot 
and thus determines 



ay X (mN/m) 


Figure 11. Effective surface tensions difference A'yLv versus 
cjyR*m (see text). The line is a linear fit for which = 0.85 
and the prefactor is 2.8. 


B. Influence of elasticity 


Results of figure clearly show that G' has a strong 
effect on the shape of the stretching-compression cycle. 
Indeed the elastoplastic force depends on the elastic de¬ 


formation of the bridge (see part IV A), and especially at 
changes of deformation direction. 

To analyse this effect in a more systematic way, the dif¬ 
ference between the force F — IV of the first step of the 
second stretching and the force corresponding to a com¬ 
pressed bridge at the same L (see figure 12 1 ) was mea¬ 
sured on each force-L plot. This force difference AF is 
plotted as a function of an estimated elastic force AFestim 
in figure [T^ . In the filament approximation this elastic 
force corresponds to: 


AFestim — 3G^ X — X 
h 

where 3G' is an estimation of the Young modulus of the 
gel, ^ is the relative variation of the bridge height h on 
the first step of the second stretching and is the 

section of the bridge at the neck after the first step of 
the second stretching. 

These two forces are nearly equal, which evidences 
the elastic behavior at the change between compression 
and stretching. This confirms that the shape of the 
stretching-compression cycle is driven by the elasticity 
modulus G'. 




Figure 12. a) Definition of AF in a force-L plot, b) Calculated 
elastic force (see text) AFestim versus AF (measured). The 
red line is a linear fit with a slope of 0.67. Inset : same plot 
in log log scale. The dashed line indicates the measurement 
error on the x-axis. On the y-axis, the error is around 25% of 
the values. 


C. Effect of the initial stress 


Many experimental force-L plots show a positive y- 
intercept for the stretching part (red fits), whereas this 
feature does not appear in the model, where this y- 
intercept is always negative. This can be explained as 
follows: the model assumes that the initial stress is max¬ 
imum (ie. the fluid has reached the yielding point) before 
the first stretching phase. But this cannot be checked ex¬ 
perimentally. It is likely that in some experiments the ini¬ 
tial stress is not maximum, although the fluid is strongly 
stretched before the beginning of the experiment. 

Three different cases are illustrated in figure which 
shows results from the model with different initial condi¬ 
tions: the full red symbols stand for a stretching phase 
beginning with a maximum stress, the pink symbols for 
the same stretching series with a zero initial stress and 
the empty symbols for a stretching phase with maximum 
negative stress (after a full compression for example). 
The ratio G'jay is set to 0.5 in this figure and the ge¬ 
ometry is filament-like. In the case of a zero initial stress 
it is clear that even if the points seem nearly aligned, a 
positive y-intercept arises and the slope decreases. 

The model shows that above a ratio G' jay « 8, the 
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Figure 13. Stretching phase in the filament geometry, with 
G'/ay = 0.5 and three different initial conditions. Full red 
symbols: maximum initial stress ao = -\/3ay. Pink symbols: 
zero initial stress ag = 0. Empty symbols: negative initial 
stress (Jo = — \/3cv (after full compression). Each set of 
points is shown with its linear fit. 


force-L curves are not sensitive any more to the initial 
stress (see figure for example) and the y-intercept of 
the stretching branch (red) is always negative. 

For the carbopol samples we used, the ratio G' jay is 
comprised between 2.3 and 6.1. Therefore the curves 
are sensitive to the initial conditions. This explains 
the occurrence of positive y-intercepts for the stretch¬ 
ing branch and also the dispersion in the limiting slopes 
values. However for our samples where G' jay > 5 (HS 
carbopols), we could obtain force-L plots with negative 
y-intercept and limited influence of initial stress, as pre¬ 
dicted by the elastoplastic model (figure [l4| . 



Figure 14. Force-L plot for carbopol 0.75% (HS) with 
G'/ay = 5.1. From the second stretching branch (empty tri¬ 
angles) we see that the initial condition is quickly forgotten. 
The mean slope is 59.5 mN/m. 


D. Surface tension 

Estimation with the bridge tensiometer. As can be 
seen in figure [^, for vanishing yield stress the upper and 
the lower slopes both tend to around 63 mN/m. This 
suggests that the surface tension of carbopol gels is close 
to this value. For higher values of the yield stress, our 
experiment clearly shows that the way an experiment 
is performed (ie. the fact that the yield stress fluid is 
stretched or compressed) influences a lot the value of the 
surface tension found via this experiment. 

The bridge tensiometer setup provides an easy way to 
measure the surface tension of yield stress fluids with 
G' ~ ScTy or more: on a force-L plot the data align on 
two limiting curves which are symmetrical with respect 
to F —W = 7 LyL. The true value of the surface tension 
is thus the mean of the slopes of the two linear fits. With 
the model, taking G' jay = 8 in the filament geometry, 
the surface tension value could be recovered in this way 
within less than 1%. Note that the condition G' jay > 8 
is usually met in a large range of yield stress fluids like 
emulsions m, clay suspensions [55] and microgels pastes 

m- 

By this method and with our HS samples for which 
G' jay > 5 we obtained as a maximal value for the mean 
slope 63.1 mN/m for 0.25% carbopol (ay = 4.6 Pa) and 
59.5 mN/m for 0.75% carbopol (ay = 15.6 Pa). 

Ascending bubble measurements. Our results with 
the bridge tensiometer are confirmed by other experi¬ 
ments with an ascending bubble setup (Teclis Tracker). 
Here again the apparent surface tension depends on the 
flow history. 

The surface tension between carbopol and air was mea¬ 
sured by injecting an air bubble in a large volume of very 
low yield stress carbopol (~ 1 Pa) and analyzing the bub¬ 
ble profile. The measurement can be static or dynamic. 
In the dynamic case, a given volume of air is injected in 
the fluid (corresponding to a given area Aq), and then a 
fixed interface area is imposed, either greater or smaller 
than Aq . The area remains then fixed thanks to a retroac¬ 
tion loop during the whole measurement, which lasts for 
about 10 minutes each time. 

The evolution of the apparent surface tension is plotted 
in figure A slow relaxation can be observed after the 
increase (resp. decrease) of the area of the bubble, at 
t = 0. The end value of the apparent surface tension is 
rather stable, but depends on the quantity of area added 
(resp. removed) at the beginning. We interpret this as 
an elastic effect which cannot relax totally because of the 
yield stress of the fluid. For a 0.25%(MS) carbopol of 
yield stress 0.3 Pa, the apparent surface tension after 10 
minutes is between 61 mN/m and 66 mN/m, depending 
on the history of the bubble. For a 0.5%(MS) carbopol 
of yield stress 1.75 Pa, the apparent surface tension is 
between 59 mN/m and 65 mN/m. 

The ascending bubble commercial device can however 
not be used for determining the surface tension of a wide 
range of carbopols, since it is not powerful enough to 
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Figure 15. Evolution of surface tension measured by an as¬ 
cending bubble method, for a 0.25%(MS) carbopol of yield 
stress ay ~ 0.3 Pa, after different changes in interface area at 
t — 0. Before t = 0 the volume of the bubble is 13 fiL and its 
area is slightly less than 25 mm^. 


push a bubble in the liquid when the yield stress is over 
2 Pa. 

Comparison with values found in the literature. Dif¬ 
ferent values can be found for carbopol surface tension in 
the literature. First, Hu et al. |29| found by a maximum 
bubble pressure method that neutralized carbopol from 
0.025% to 0.1% had the same surface tension as pure wa¬ 
ter on a wide temperature range. 10 years later, Manglik 
et al. EOl found by the same method, but using dynamic 
and static measurements, that equilibrium surface ten¬ 
sion of carbopol decreased clearly when concentration 
increased. For their maximum polymer concentration 
(0.2%), the surface tension is measured at 69 mN/m. 
However their carbopols seem not to be neutralized and 
the rheological data are unclear. An explanation for the 
discrepancy between these references’ results and ours 
could be that the yield stress and elasticity effects were 
not taken into account, even if Manglik et al. corrected 
their dynamic measurements with a viscosity term. 

More recently, Boujlel and Coussot [5] used a plate 
withdrawal method to study the effect of the capillary 
Bingham number = cry x D/jj^y on the withdrawal 
force, changing the yield stress and the dimension D of 
the plate. An extrapolation to Be = 0 allowed them to 
estimate the surface tension of carbopol to 66 mN/m. 


Geraud et al. E] performed capillary rises of neutral¬ 
ized carbopol gels and could extract a value for the sur¬ 
face tension from the maximum rise height as a function 
of the gap. They found cos0o = 48±3 mN/m, which 
cannot be compatible with Boujlel’s results unless the 
contact angle Oq is at least 43°. Yet the contact angle 
of carbopol droplets freshly dropped on the same sur¬ 
faces used by Geraud et al. was always measured smaller 
than 25°. This value would imply a maximal value of 53 
niN/m for the surface tension of carbopol. 

Our work is able to reconcile these different measure¬ 
ments, as it shows that the yield stress and elasticity 
of a fluid give rise to a separation in two different val¬ 
ues for the effective surface tension, depending on the 
flow history. The maximum bubble pressure method and 
the plate withdrawal method being always performed in 
surface extension, the energy cost comes from surface in¬ 
crease as well as from the rheological resistance of the 
gel. We thus expect a higher value than for capillary 
rises where the fluid is moving forward driven by surface 
tension. 


VI. CONCLUSIONS AND PERSPECTIVES 

In this work we have shown that the measurement of a 
yield stress fluid surface tension is strongly influenced by 
the fluid rheology and the protocol. In particular, the di¬ 
rection of the flow during or even before the measurement 
gives rise to excess forces, either positive or negative, that 
are difficult to compute in a general case. It is thus useful 
to perform measurements with a setup allowing to test 
different flow histories. 

We have proposed a method to measure the surface 
tension of yield stress fluids of high elastic modulus (G' > 
8 cry). The surface tension value is given by the mean of 
the two limiting slopes in a force-L plot. 

To go further, capillary bridges are a convenient tool 
to study the surface tension and the rheology of simple 
or complex fluids. Our method can be very precise when 
it comes to measure the surface tension of simple fluids. 
With complex fluids, the complexity of the real shape 
of the bridge does not change qualitatively the adhesion 
force compared to model geometries. 

Future work will focus on the influence of wall slip and 
line pinning on the adhesion force of capillary bridges 
and more generally on the measurement of surface ten¬ 
sion. In this scope, other experiments, including surface 
fluctuation spectra analysis m, have already started. 
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APPENDIX A: TYPICAL FORCE EVOLUTION 
DURING AN EXPERIMENT 


The evaporation of water and carbopol is not negligible 
during a whole experiment (about 30 minutes). This is 
reflected in the force evolution in time even in the absence 
of flow (figure 161. The force decreases slowly, at an 
approximative rate of 1 fjiN per second. 

To minimize the uncertainties on the force, its value 
is saved only when this slow evolution is reached and 
the picture is taken at the same moment. The resulting 
uncertainty is very small compared to the force values 
(100 to 1000 /xN) and it is taken into account when fitting 
the force-L plots. 



Figure 16. Force evolution in time during a typical stretching 
step. The red ellipses represent the moments when the force 
value and the picture are saved. 


APPENDIX B: CALCULATION DETAILS IN THE 
PANCAKE GEOMETRY 


The differences between the filament and the pancake 
geometries lie in the expression of the geometrical pa¬ 
rameter L, the dominating terms of the stress tensor and 
then the calculation of the normal force. 

We recall the approximate expression of L in a quasi- 
cylindrical geometry where the volume is U = 


L « 



2V cos 00 


In the case of a very fiat drop, the dominant curvature 
is in the (r, z) plane so that L rv 2V cosOo/h^- 

In addition the deformations and dissipation are dom¬ 
inantly due to shear along the z direction, in the lubri¬ 
cation approximation. In this case, we cannot use a uni¬ 
form description like in the filament geometry, but we 
need to describe the stress profile at the wall across the 
axis i7wan(?') = = Zwaii)- Note that the region 

of high stress in the pancake geometry is located at the 
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walls whereas the highest stress region is located in the 
neck of the filament. 

The extra pressnre due to the elastoplastic flow is here 
denoted p = Ap — ApLapiace- The stress balance reads: 

darz dp 
dz dr 


We also re-evaluate the values of the radial positions r' = 

r /\h 

r -and finally evaluate at each step the traction 


force 


2h 


cRn 


Fe-p = 


2 2^wall(^) , 
7rr -;- dr 


which is a function of r only, so that — = —;— . This 

dr h 

allows for the determination of the elastoplastic traction 
force: 


p — 
ep — 


Each step starts with a given stress profile cro(r). For a 
variation in height Ah, the corresponding step in elastic 

T 

stress is Aa(r) = SG'Ah—. For each position, we thus 
evaluate the new stress value at the wall: 

{ —ay if o'o(r) -f Acr(r) < —ay 

cro(f) + Acr(r) if-cry < cro(r) + Acr(r) < cry 
-l-cry if cro(r) -I-Acr(r) > -l-cry 


oRn 


27rrp{r)dr 


2dp, 
Trr -—dr 
dr 


2 wall 7 

Trr - dr 

h 


In practice the drop is assumed to be initially stretched 
just at the onset of yielding so the initial stress is set at 
cro(r) = ay(r / Rjq) (this linear variation corresponds to a 
pure elastic deformation). Finally we stretch then com¬ 
press and finally stretch again the drop from h = 0.5 mm 
to h = 1.5 mm by successive steps Ah = 0.1 mm. For 
each step, the total traction force is the sum of the capil¬ 
lary force and the elastic one. Typical curves of traction 
forces as a function of L are plotted in figure The 
same phenomenology is observed in the pancake and fil¬ 
ament geometries: stretching (resp. compression) of the 
capillary bridge is associated with an increase (resp. de¬ 
crease) of the apparent surface tension, which is more 
pronounced for higher yield stresses. The force cycles also 
depend on the elastic modulus and initial stress condi¬ 
tion. A noticeable difference between the two geometries 
is the explored range of the geometrical parameter L: the 
pancake geometry corresponds to more compressed drops 
characterised by higher L values than in the filament ge¬ 
ometry. 







